Modelling the Nonlinear High-Prequency Response of a Short Josephson Junction 

under Two-Frequency Irradiation 
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The nonlinear response of a short Josephson Junction (JJ), 
being irradiated simultaneously with two high-frequency sig- 
nals, has been studied in the framework of the nonlinear 
Resistively-Shunted Junction (RSJ) Model. One of the sig- 
nals, hereafter referred to as "probe signal" , has a small am- 
plitude Ip r < I c (J c is the critical current of the JJ) and 
frequency f pr , and is used to monitor the response of the 
junction to the other high-power signal with amplitude I prn 
and frequency f pm , hereafter referred to as "pump signal". 
Varying the frequency ratio f pm / fpr from 0.5 to 100, and 
the current amplitude of the probe signal from 0.01 to 0.9 of 
I c , we found that the dependence of the junction impedance 
at the frequency f pr , Zg, versus I pm preserves its general 
features, independent of fpm./ fpr and I pr /Ic values. At the 
same time, some particular features, like negative values of 



Re(Zs) and "fine" structure of the steps in Zi pT (I pm ) are 
observed for fpm,/ fpr < 1 and for particular values of I pr /I c . 
In general, the behavior of Z a pr (I P m) is rather different from 
that predicted by the nonlinear RSJ model for a short JJ in 
the regime of single- frequency irradiation, when one and the 
same signal plays the roles of the pump and the probe signals 
simultaneously. Possible applications of the model are briefly 
discussed. 
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I. INTRODUCTION 

The high-frequency nonlinear response of a Josephson 
Junction (JJ) is of significant interest because JJs are nec- 
essary units of almost all active microwave and millimeter 
(mm) wave devicesu, and because weak links, which are likely 
to be present even in the highest quality samples of high-T c 
superconductors (HTS), can be modelled as JJsB. The RSJ 
model it. often used to simulate the characteristics of JJ-based 
devicesa, and was shown to give a good agreement with ex- 
perimental data on point-contacts and microbridges in the 
microwave and mm wave ranges, where the capacitance of 
the junction can be neglected. 

In the present paper we report simulation of the surface 
impedance of a JJ at the frequency of the low amplitude 
signal, hereafter referred to as "high-frequency response", as 
a function of the current amplitude of the other elevated- 
power high-frequency signal. This case can be considered as 
a model of a microwave-biased electromagnetic radiation de- 
tector, which has been shown to have an improved sensitive 
and noise figures when compared with dc-biased detectorsu 
The other possible implication of the model is modeling the 
performance of microwave parametric amplifiers and Joseph- 
son mixerail. In addition, the above model can also describe 
the nonlinear microwave response of superconducting weak 



links, which is offamnriinvestigated using the so-called "pump- 
probe" techniqueo'S'u. This technique is of a particular inter- 
est because it allows one to modulate or to pulse the powerful 
microwave signal, whilst measuring the surface impedance of 
the sample with the help of the other low-amplitude contin- 
uous wave microwave signal at a different frequency. In such 
a way, the pump-probe method avoids substantial heating ef- 
fects and allows the study of intrinsic nonlinear phenomena 
in superconductors. However, using the pump-probe tech- 
nique, one has to know how the nonlinear surface impedance, 
measured at the pump frequency, relates to that measured at 
the frequency of the probe signal, with respect to which the 
superconductor is in the linear regime. Although the model 
proposed in this paper cannot be considered as a model of 
the nonlinear response of HTS, an extension of the model 
to the case of a 2D JJ array with a random distribution of 
Ic • -Rn-products (as cecently proposed by Herd et al. for the 
single-frequency caseO) would allow the description ftfpfl^per- 
iments on HTS thin films in the pump-probe regime 



II. NUMERICAL SIMULATION 

In the case of two signals applied to a short JJ, the sine- 
Gordon Equation (which describes time dependence of the 
order parameter phase difference ip across the junction) in 
dimensionless form can be written as followsEj 



dr 



, sin(fi pm r) + ipr sin(f2 pr r) — sin ip, 



(1) 



where r = fit, i pm = I pm /Ic 



— Ipr/Ic 



Ipm and Ipr 



current amplitudes of the pump and probe signals respec- 
tively, Upm = 2nf P m. and Lu pr = 2nf pr the corresponding 
circular frequencies, I c the critical current of the junction, 



fi = 2eR n I c /h, fi p 



r /f3, and R n surface 



resistance of the junction in the normal state. Equation (|lj) 
does not take into account the effect of the junction capaci- 
tance, which can be neglected at microwave frequencies. 

Generally, the solution of (jl|) is not necessarily periodic 
with time. Only in the case when the ratio Q pm /^pr is an 
integer is the solution of ([!]) periodical with a period equal to 
the least common multiple (LCM) of the pump and the probe 
signal periods. If we expand the phase derivative <p into a 
Fourier series with respect to time over the LMC period of the 
two signals, we obtain coefficients which couple the voltage 
~ <p to the current I r f. If then we single out the series' terms 
at the probe frequency, we obtain the surface impedance at 
the relevant frequency as follows 



Z s = R s + jX s = lim 



Q,Rn 



<p(r) exp(jfl pr .r)dr, (2) 



1 




FIG. 1. Normalised surface resistance R a /R„ (solid line) 
and surface reactance X„/R n (broken line) for a short Joseph- 
son junction, as a function of normalised pump current am- 
plitude ipm = Ipm/Ic, simulated within the framework of the 
RSJ model for the two-frequency case. Here, the frequency 
ratio SI 

each of the figures 



pm/Qpr = 2, and the probe current i pr is that given in 



where Q = 2tt/T, and T is the LCM of the pump and probe 
signal periods. Because one can assume different initial con- 
ditions for (|l]), its solution is not strictly periodic with T, and 
hence in (Q) integration over a few periods is required to get 
an appropriate convergency of the results. 

The parameters used for simulation are as follows: R n = 



l(T d «, Ic 
7.2)-10 10 Hz, 



0.5 A, fpn 



3.6 • 10 iu Hz, f P r =(0.036- 



n prn = 2TYhfp m /(2eI c R„) = 0.141. 



pm /S^pr 



Results of the simulation for different ratios of Q, 
varied from 0.5 to 100, and for different values of the probe 
current amplitude Lp r (0.05, 0.25, and 0.7) are plotted in 



Fig. 
Fig. 



Fig. |, Fig^ § and Fig. | For Q pm /n pr > 1 (Fig. Qa, 
2a and Fig. |3ja) at low i pr (0.05), the surface resistance 



FIG. 2. Rs/Rn and X s /Rn for a short Josephson junction, 
as a function of i pm = Ipm/Ic, simulated within the frame- 
work of the RSJ model for the two-frequency case. Here, the 
frequency ratio fl pm /0,p r = 10, and the probe current i pr is 
that given in each of the figures. 



rather different from that expected for the single-frequency 
situation. In the latter case, the surface resistance R s has a 
staircase-like form, starting to increase rapidly for i pm > 1, 
and gradually approaching the R n value With 
regards to the surface reactance X a , for the single-frequency 
case it oscillates around ^ero with amplitude decaying with 



increased i pm (see, e.g. 



However, in the two- frequency 



regime, no obvious decay of the peaks amplitude in X a up 
to i pm — 4 is observed. In addition, every oscillation peak 
seen in X s (i V m) in the single- frequency regime translates into 
a peak with a complicated structure, containing many up- 
ward and downward minor peaks of smaller amplitude. With 
increased i pr the major peaks in xl pT (i pm ) broaden, and 
more complicated "fine" structure of minor peaks develops 
(see (Fig. l]b,c, Fig. ^b,c and Fig. ^|d,c)) . As far as R. 



Rl pr as a function of I pm is a sequence of peaks of about concerned, an increase of i pr leads to the appearance of steps 
the same height with a "pedestal" level approximately equal 
to the value of R{ pr in the linear regime. This picture is 



in R a pr (i P m), similar to those seen in R s (i P m) for the single- 
frequency regime. The higher i pr , the more steps are observed 
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FIG. 3. Rs/Rn and X s / R n for a short Josephson junction, 
as a function of i pm = Ipm/Ic, simulated within the frame- 
work of the RSJ model for the two-frequency case. Here, the 
frequency ratio il pm /Q pr = 100, and the probe current i pr is 
that given in each of the figures. 
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FIG. 4. Rs/Rn and X a /Rn for a short Josephson junction, 
as a function of i pm = Ipm/Ic, simulated within the frame- 
work of the RSJ model for the two-frequency case. Here, the 
frequency ratio Q pm /0, pr = 0.5, and the probe current i pr is 
that given in each of the figures. 



in Rs pr (i pm ), before it levels off and starts to oscillate near 
some average value around unity for i pm S> 1. 

Contrary to the case of fip m /f2 pr > 1, when the major 
peaks in R B pr {i P m) are almost symmetrical with respect to 
a vertical line drawn through the middle of their width, in 
the case of Sl„ m /Sl pr < 1 these peaks are obviously asymmet- 
ric (see Fig. kf). Another distinctive feature for this case are 
discontinuous double-peak structures (upward peak followed 
by downward peak) at the beginning and the end of every 
major peak. One further significant difference of R a pT (i pm ) 
in the case of Q. pm /£l pr < 1, as compared to the case with 
f2 pm /fi pr > 1, is the appearance of regions with negative val- 
ues of R s pr . This means that under these particular condi- 
tions the JJ can contribute energy to the external circuit, i. e. 
it works as a generator. This effect was theoretically predicted 
and experimentally observed in JJs made of low-temperature 
superconductors, and was called "the effect-pf nondegenerated 
single-frequency parametric regeneration"!!! As the theoret- 



ical analysis showed, this phenomena can be realised in any 
parametric element, a reactive parameter of which can take 
negative values with changing timaij. 

All other features of R{ pT (i pm ) for the case of Q pm /Q pr < 1, 
such as the appearance of steps, an increase in their number, 
and a shift of the oscillatory part of the dependence to higher 
i pm with increased i pr , are similar to those seen in the case of 
Fl P m/Q pr > 1. As far as X s pT (i pm ) is concerned, features like 
asymmetry of the major peaks and discontinuous double-peak 
structures are observed, similar to those present in R s p r (i pm )- 

A. Implications for applications 

One of the possible applications of the two-frequency 
regime, simulated in this paper, is a microwave-biased JJ de- 
tector. An advantage of this regime is that the amplitude of 
the oscillation peaks (or, equivalently, the impedance steps) 
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can be made several times (up to a factor of 3-4) higher than 
the resistance of the junction in the normal-state, especially 
at a small probe current (see Fig. Ek, Fig. ^a and Fig. ^a). 
This should lead to enhanced sensitivity of the detector, as 
compared to the single-frequency regime. Despite that sim- 
ilar results (increased step heights) have been also obtained 
for a dc-biased JJu, a microwave-biased detector was shown 
to benefit from reduced noise temperature andi-enhanced re- 
sponsivity as compared to the dc-biased oneau. Moreover, 
since the step amplitudes in the saturation regime (I pm 2> 1) 
are almost independent of the pump current, such a detector 
will possess amplitude-independent sensitivity. 

In the case when Q pm /£l pr < 1, the J J can be used as a 
parametric amplifier operating in the single-frequency nonde- 
generate regime which, when compared to the self-pumped 
regime, was shown to give a reduced noise temperature and 
a narrows* frequency response of the junction at the probe 
frequency!], allowing one to use high quality factor resonators 
for matching the junction with the external circuits. 

III. CONCLUSION 

The numerical simulation performed by us has shown that 
the nonlinear high-frequency response z{ pr of a short JJ in 
the regime of two-frequency irradiation can be rather differ- 
ent from the surface impedance Z 3 measured in the single- 
frequency regime. Depending on the ratio of the pump to 
the probe frequencies, a number of new features in Z s pr (i pm ) 
are predicted. Among them are the absence of the steps 
in R{ pr '{i P m) at low probe currents (i pr < 0.05); persistent 
oscillations of R s pr (i P m) around some average value which 
tends to unity with increased i pm ; a multiple-peak structure of 
X B pr (ipm), which becomes more complicated with increased 
ratio Q P m/Qpr; appearance of regions with negative values of 
surface resistance in R a pr (i pm ) for the case of Q, prn /fl pr < 1. 

At the same time, there are some features which are similar 
to the those in the single-frequency regime. These are the 
appearance of steps in R{ pl ~ {i pm ) with increased i pr , and the 
oscillation of X s pr (i pm ) around a zero- Xi pr value. 

The model presented here was shown to give a useful basis 
knowledge for an application of the JJ as a microwave- biased 
detector of electromagnetic radiation. 
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